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ABSTRACT 

We use gravitational microlensing of the four images of the z = 0.658 quasar RXJ 1 131-1231 to measure the 
sizes of the optical and X-ray emission regions of the quasar. The (face-on) scale length of the optical disk 
at rest frame 400nm is R\,o = 1.3 x 10 15 cm, while the half-light radius of the rest frame 0.3-17 keV X-ray 
emission is R\/2,x = 2.3 x 10 14 cm. The formal uncertainties are factors of 1.6 and 2.0, respectively. With the 
exception of the lower limit on the X-ray size, the results are very stable against any changes in the priors used 
in the analysis. Based on the H/3 line-width, we estimate that the black hole mass is Mini ~ 10 8 M Q , which 
corresponds to a gravitational radius of r g ~ 2 x 10 13 cm. Thus, the X-ray emission is emerging on scales of 
~ 10r g and the 400 nm emission on scales of ~ 70r ? . A standard thin disk of this size should be significantly 
brighter than observed. Possible solutions are to have a flatter temperature profile or to scatter a large fraction 
of the optical flux on larger scales after it is emitted. While our calculations were not optimized to constrain the 
dark matter fraction in the lens galaxy, dark matter dominated models are favored. With well-sampled optical 
and X-ray light curves over a broad range of frequencies there will be no difficulty in extending our analysis to 
completely map the structure of the accretion disk as a function of wavelength. 

Subject headings: accretion - accretion disks - black hole physics - gravitational lensing — quasars: individual 
(RXJ 1131-1231) 



1. INTRODUCTION 

A significant problem for theoretical studies of quasars is 
that we cannot spatially resolve their emission regions to test 
models (e.g. Blaes 2004). For example, in this paper we study 
the gravitational lens RXJ 1231-1131 (RXJ1131 hereafter), 
where we observe four images of a z s = 0.658 quasar lensed by 
a zi = 0.295 elliptical galaxy (Sluse et al. 2003). Based on the 
H/3 line-width from Sluse et al. (2003), and a magnification 
corrected estimate of the continuum luminosity, we estimate 7 
that the black hole mass, M BH , is M\ m = ( 1 .3 ± 0.3) x 10 8 M Q . 
This corresponds to a gravitational radius of 
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(1) 



that subtends only 0.001/T 1 micro-arcseconds. 

Gravity, however, has provided us with a natural telescope 
for studying the structure of quasars through the microlens- 
ing produced by stars in the lens galaxy (see the review by 
Wambsganss 2006). Microlensing has a natural outer length 
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scale corresponding to the Einstein radius of the stars, 
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where (M) is the mean mass of the stars, and the distances 
Dql, D os and D LS are the angular diameter distances between 
the observer, lens and source. The microlenses also gener- 
ate caustic lines on which the magnification diverges, which 
means that our gravitational telescope can, for all practical 
purposes, resolve arbitrarily small sources. The size of the 
source is encoded in the amplitude of the microlensing vari- 
ability as the source, lens, and observer move relative to the 
caustic patterns - big sources have smaller variability ampli- 
tudes than small sources. The technique can be applied to any 
emission arising from scales more compact than a few (Re)- 

If we model the accretion disk by a thermally radiating thin 
disk with a temperature profile of T oc R~^I A (Shakura & Sun- 
yaev 1973) 8 , we can measure the scale R\ defined by the 
point where the photon energy equals the disk temperature, 
kT = hc/Xrest, by two routes other than microlensing. First, 
we can estimate it from the observed flux at some wavelength. 
For example, at I-band the radius is 
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where r# = c/Hq is the Hubble radius, and i is the in- 
clination angle of the disk. Based on HST observations 
(Sluse et al. 2006; Kozlowski et al. 2009), we esti- 
mate that the magnification-corrected flux is I ~ 20.7 ±0.1 
mag (X bs = 0.814/j,m), which corresponds to an R-band 
(400 nm in the quasar rest frame) size of Rx.o = (3.5 ± 0.2) x 

8 In our present analysis we can neglect the drop in temperature and emis- 
sion near the inner edge of the accretion disk as it has little effect on the 
results. 
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10 14 (cos/) x l 2 cm or about 18^. The flux size depends on 

the mean magnification of the images as 1/ J (/j), which can 
introduce a ~ 50% systematic uncertainty into this size esti- 
mate. Second, thin disk theory predicts that 
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which implies an R-band disk size R\ o = 1 .6 x 10 15 cm (82r ? ) 
if the disk is radiating at the Eddington limit {L/L E ) = 1 with 
an efficiency of 77 = 0. 1 . Note that these two size estimates can 
be reconciled only if (L/rjL E )(M BH /M m \) 2 — 0.1(cos/)~ 3//2 , 
corresponding to a sub-Eddington accretion rate, an overesti- 
mated black hole mass, or a problem in the disk model since 
there is no evidence for the 1-2 mags of extinction in the lens 
galaxy that would be needed raise the flux size up to that from 
thin disk theory (Eqn.|5]l. Adding the inner disk edge or us- 
ing a simple relativistic disk model (Novikov & Thorne 1973, 
Page & Thorne 1974) changes this problem little. 

The expected size of the X-ray emitting regions is more 
problematic because there is no comparably simple model 
for our theoretical expectations. There is a general consen- 
sus that the X-ray continuum emission is due to unsaturated 
inverse Compton scattering of soft photons by hot electrons in 
a corona surrounding the inner parts of the accretion disk (see 
the review by Reynolds & Nowak 2003), but the extent and 
geometrical configuration of the X-ray emission region is an 
open question. The X-ray continuum from the corona illumi- 
nates the disk to produce Fe Ka emission lines, whose broad 
widths indicate that they are generated close to the inner edge 
of the accretion disk (e.g. Fabian et al. 2005). 

While there were a number of early attempts at estimating 
accretion disk sizes using microlensing (e.g. Wambsganss, 
Schneider & Paczyhski 1990, Rauch & Blandford 1991, 
Wyithe et al. 2000b, Wambsganss et al. 2000, Goicoechea 
et al. 2003), it is only in the last few years that it has be- 
come possible to make large numbers of microlensing size 
estimates. In particular, Pooley et al. (2007) argue that the 
optical sizes estimated from microlensing must be consider- 
ably larger than the optical "flux" sizes of Eqn. [3] This was 
confirmed by Morgan et al. (2009) in a more detailed analy- 
sis that also found that the optical sizes agree better with the 
thin disk size estimate (Eqn. |5]i than the flux size and have a 

scaling with black hole mass consistent with the M 2 J^ scaling 
for Eddington-limited thin disks. 

Recent studies have started to examine the temperature de- 
pendence of disks through the scaling of disk size with wave- 
length (Anguita et al. 2008, Poindexter et al. 2008, Agol et 
al. 2009, Bate et al. 2009, Floyd et al. 2009, Mosquera et 
al. 2009). Studies of the microlensing of the X-ray emission 
are more limited, but indicate that the X-ray emission is much 
more compact than the optical (Dai et al. 2003, Pooley et al. 
2006, 2007, Kochanek et al. 2006, Morgan et al. 2008, Char- 
tas et al. 2009), tracking much closer to the inner edge of the 
accretion disk. In this paper we estimate the sizes of the opti- 
cal and X-ray emission regions of RXJ1 131 using microlens- 
ing. In §2 we describe the data and the analysis method. In 
§3 we discuss the results, their implications and directions for 
further research. We use an Q,q = 0.3 flat cosmological model 
with H = lOO/i km s" 1 Mpc" 1 and h = 0.7. 



2. DATA AND ANALYSIS 

The optical data consist of the five seasons of R-band mon- 
itoring data described in Kozlowski et al. (2009). For our 
present analysis we simply shifted the light curves by their 
measured time delays (Kozlowski et al. 2009). The X-ray 
data, all ACIS observations from the Chandra Observatory, 
consist of the epoch presented by Blackburne et al. (2006) 
plus the 5 epochs presented in Chartas et al. (2009). Each of 
the Chartas et al. (2009) epochs consisted of a 5 ksec obser- 
vation using ACIS-S3 in 1/8 sub-array mode from which we 
measure the 0.2-10 keV flux. Chartas et al. (2009) also reana- 
lyzed the Blackburne et al. (2006) data to properly correct for 
the "pile-up" effect. We do not use the X-ray fluxes of image 
D in our analysis because we cannot presently be certain its 
flux ratios relative to A-C are unaffected by source variabil- 
ity given the roughly 3 month time delay between D and A-C 
(Kozlowski et al. 2009). As we can see from Fig. [5] the X-ray 
source must be more compact than the optical source because 
the X-ray flux ratios are dramatically more variable. 

A full description of our microlensing analysis method is 
presented in Kochanek (2004) and Kochanek et al. (2006). 
In essence, we create the microlensing magnification patterns 
we would see for a broad range of lens models and source 
sizes, then randomly generate light curves to find ones that fit 
the data well. We then use Bayes' theorem to combine the re- 
sults for the individual trials to infer probability distributions 
for physically interesting variables including the uncertainties 
created by all the other variables. 

We fit the lens as in Kozlowski et al. (2009), modeling it as 
a R e = 1 '.' 7 de Vaucouleurs model for the stellar distribution 
embedded in an NFW halo. We consider a sequence of mod- 
els described by /», the fraction of mass in the stellar com- 
ponent relative to a constant mass-to-light ratio model with 
/„ = 1 and no halo. We include models with /» = 0.1 to 1 in 
equal steps, and the time delay measurements favor /* ~ 0.2. 
These lead to the values for the convergence k, shear 7 and 
fraction of the convergence in stars k*/k reported in Table. 1. 

The stars creating the microlensing magnification were 
drawn from a power law mass function dN/dM oc M~ 13 with 
a ratio of 50 between the minimum and maximum masses that 
roughly matches the Galactic disk mass function of Gould 
(2000). We know from previous theoretical studies that the 
choice of the mass function will have little effect on our con- 
clusions given the other sources of uncertainty (e.g. Paczyhski 
1986, Wyithe et al. 2000a). The mean mass (M) is left as 
a variable with a uniform prior over the mass range 0.1 < 
(M/M Q ) < 1.0. 

For each model we generated 8 random realizations of the 
star fields near each image. The magnification patterns had an 
outer scale of 10(R E ) = 4.6(M/M ) 1/2 x 10 17 cm and a pixel 
scale of 10(/? £ )/8192 = 5.6(M/M ) I/2 x 10 13 cm ~ 3r g , so 
we should be able to model sources as compact as the inner 
edge of the accretion disk. We modeled the relative velocities 
as in Kochanek (2004), where for RXJ1 131 the projection of 
the CMB dipole velocity (Kogut et al. 1993) on the lens plane 
is 47 km/s, the lens velocity dispersion estimated from the 
Einstein radius is 350 km/s, and the estimated rms peculiar 
velocities of the lens and source galaxies are 180 and 140 km/s 
respectively. 

The source model for both the optical and X-ray sources is 
a face-on disk with a temperature profile T oc i?~ 3/4 radiating 
as a black body (Shakura & Sunyaev 1973), so the surface 
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brightness profile of the disk is 
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with the single parameter being the scale length R\. While 
it is true that this profile lacks the central drop in emissivity 
and that it is not a physical model for the non-thermal X-ray 
emission, the microlensing analysis is not sensitive to these 
details. The estimate of the half-light radius (R\/2 — 2A4-R\) 
is essentially independent of the assumed profile (Kochanek 
2004, Mortonsen, Schechter & Wambsganss 2005). We used 
a 46 x 61 logarithmic grid of trial source sizes for the X-ray 
and optical sources with a spacing of 0.05 dex. 

We do, however, allow for the possibility that fraction 
fnoii = to 40% of the optical emission is generated on scales 
much larger than the disk and is unaffected by microlensing. 
Such large scale emission could have two physical origins. 
First, the optical continuum can be significantly contaminated 
by emission lines, both the obvious broad lines and the less 
obvious Fe and Balmer pseudo-continuum emission (~ 30% 
of the emission in some Seyferts, Maoz et al. 1993), that are 
believed to be produced on much larger scales than the disk. 
For our R-band light curves, there are no strong emission lines 
in the filter band pass, but the blue edge of the Balmer con- 
tinuum emission (~ 6000A) does lie inside the band pass 
(roughly 5700-7200A). Second, even if the observed pho- 
tons were generated by the accretion disk, a fraction could 
be scattered on much larger scales, leading to an effectively 
larger source. These two possibilities are not equivalent, as 
the line emission is due to reprocessing of shorter wavelength 
UV photons rather than the observed R-band continuum. 

A basic problem for any microlensing analysis is the degree 
to which the "macro" lens model correctly sets the average 
magnifications. Each light curve, m,(f) = s(t)+ fii + 6fii(t) + A,- 
is defined by the source light curve s(t ), the macro model mag- 
nification fij, the microlensing magnification <5/i,(f) and a pos- 
sible offset Aj. These offsets can be non-zero due to problems 
in the macro model or the presence of unrecognized substruc- 
tures that perturb the magnifications (e.g. Kochanek & Dalai 
2004), because of differential absorption due to dust or gas 
in the lens galaxy (e.g. Falco et al. 1999, Dai & Kochanek 
2009), or due to contamination of the light curves by flux from 
the quasar host or lens galaxy. For the latter two possibilities, 
the offsets would differ between the optical and X-ray light 
curves. Given a sufficiently long light curve, the offsets can 
be determined from the data, but they are poorly constrained 
until the light curve is a good statistical sampling of the mag- 
nification pattern. We will consider four treatments of this 
problem to ensure that such systematic problems do not af- 
fect our results. The basic division we will refer to as Cases 
I and II. In Case I we allow the magnification offsets A,- to 
float independently for the two bands constrained by a term 
A 2 /2cr 2 in the log likelihood with a = 0.5 mag. In Case II we 
allow them to float, but use the same offsets A,- for both the 
optical and X-ray light curves. These are weak constraints, so 
the resulting distributions for the offsets are broad. To make 
sure we are not allowing too much freedom, we also examined 
limiting the range of the offsets to |A,-| < 0.3 mag in Cases I' 
and IF. 

The advantage of the less constrained strategies is that they 
are robust against the systematic errors that can plague the 
absolute magnifications of the images. It is certainly true that 
analyses using only the "DC" flux ratios (e.g. Pooley et al. 
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FIG. 1. — The probability distributions for the size of the X-ray (top) and 
R-band (400 nm in the rest frame) optical (bottom) emission regions for the 
log (solid) and linear (dashed) size priors. These sizes are marginalized over 
/nofi- The vertical lines mark the gravitational radius r g for a M1131 black 
hole, the Einstein radius for (M) = Mqand the accretion disk size estimates 
based on either the I-band flux (Eqn. [3) or thin disk theory (Eqn. |5J- The 
microlensing sizes and the I-band flux estimates can also be scaled by a 
(cos iy l I 2 inclination dependence from the assumed face-on case ((' = 0°). 

2006, 2007, Bate et al. 2009, Floyd et al. 2009) require 
less data than our "AC" approach, but they can also lead to 
conclusions dominated by these systematic errors. The "AC" 
approach also has the advantage that including the effects of 
the velocities allows us to estimate source sizes in centimeters 
without simply assuming a mean mass (M). However, when 
we use loose priors on the DC flux ratios, we lose signifi- 
cant information on the locations of the images relative to the 
magnified and demagnified regions of the patterns. As such, 
it is a conservative approach. We consider all four offset treat- 
ments in order to explore their consequences on estimates of 
the source size and the amount of dark matter in the lens. 

We used 8 statistical realizations of the microlensing mag- 
nification patterns for each of the 10 stellar surface densities 
(/*) and for 5 un-microlensed fractions of optical light (faOfi)- 
We modeled the data sequentially, making 10 6 trials for each 
optical source size and case, and then fitting each trial that 
was a reasonable statistical fit to the optical data to the X-ray 
data for each of the X-ray source sizes. In this second step we 
considered both the case where the X-ray and optical share 
the same intrinsic flux ratios and where they are allowed to 
differ. 

3. RESULTS AND DISCUSSION 

Fig. [T] shows the main result for the estimated size of the 
X-ray and optical emission regions. These combine all four 
treatments of the magnification offsets. Also note that in 
order to preserve the meaning of size ratios in Fig. Q] we 
used the scale R\ of a face-on disk for both. Physically, 
the X-ray emission is better characterized by its half light ra- 
dius, Rx/2 = 2.44R\. The scale length of the thin disk also 

scales as cos" 1 ' 2 ; if not viewed face on. We show the re- 
sults for two different priors on the disk sizes, a logarithmic 
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FIG. 2. — The correlated probability distributions for the size of the optical 
and X-ray source sizes in Einstein units of (M/Mq ) l ' 2 cm. The contours are 
drawn at the 68%, 90% and 95% maximumum likelihood contours for one 
variable for log (solid) and linear (dashed) size priors. 
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FIG. 3. — Source size dependence on parameters. The optical (top) and 
X-ray (bottom) source sizes (R\) as a function of the fraction /no/j of the 
optical flux that is not microlensed. The triangles, squares, pentagons and cir- 
cles show the results for the Case I (independent magnitude offsets for both 
bands), Case II (common magnitude offsets), Case I' (independent offsets 
limited to |A, | < 0.3) and Case IF (common offsets limited to A,- < 0.3) 
treatments of the magnitude offsets. The filled (open) symbols show the re- 
sults for the logarithmic (linear) priors on the source sizes. The horizontal 
lines show the same physical scales as in Fig.fJJand the dashed curve shows 
the expected scaling of the optical size with fno^ if we keep the half-light 
radius of the optical emission fixed. The half light radius of the disk emission 
is always R\/2 = 2.44i? A , but the half light radius of the disk emission com- 
bined with the unmicrolensed large scale emission grows with /nop, reaching 
R 



1/2 = 5.87R A for/n O/ ,=0.4. 
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FIG. 4. — Dependence on halo structure. The solid/squares and dot- 
ted/triangles show the likelihood of /» , the fraction of mass in the stellar 
component in the lens model compared to a constant M/L model (/» = 1), 
for weakly constrained (Case I+II) or strongly constrained (Case I'+IF) treat- 
ments of the magnification offsets. Dark matter dominated models are always 
favored, but the low /» models implied by the time delays are only strongly 
favored when we force the offsets to be small. The line without points shows 
the fraction 1 - k* /k of the local surface density near image A that is com- 
prised of smoothly distributed dark matter. 

(P(R\) oc 1 /R\) and a uniform (P{R\) oc constant) prior, and 
this has minor effects for the optical estimate and significant 
effects for the X-ray estimate. For the logarithmic prior we 
formally find that the (face on) optical disk scale length is 
log(fl A . /cm) = 15.11 (14.89 < log(fl A . /cm) < 15.32) and 
that the X-ray half-light radius is log(R l / 2 .x/ cm ) = 14.36 
(14.04 < logCR^x/cm) < 14.68). These estimates use both 
the prior on the velocities and a uniform prior for the mass 
over the range 0.1 < (M/M®) < 1. We will focus on re- 
sults including this mass prior, but note that if we make no 
assumption about (M), the sizes change little. With only 
the velocity priors we find \og{R\ o/cm) = 15.02 (14.75 < 
log(/? A . /cm) < 15.27) and log(fl, Az /cm) = 14.02 (13.67 < 
\og(R\/2.x/ cm ) < 14.38). The source sizes become a little bit 
smaller, but the net effect is very modest for the reason out- 
lined in Kochanek (2004). 9 

The X-ray size is more sensitive to the priors because the 
convergence of the probability distributions for small sources 
is poor when the light curve is sparsely sampled. Fig.|2]shows 
likelihoods for the source size in the Einstein units used for the 
basic calculations, and we see that they converge for small X- 
ray sizes when we use a linear prior but not for a logarithmic 
prior. The problem is not due to the pixel scale of the maps, 
but due to the lack of a well-sampled peak in the X-ray data. 
Very small source sizes are constrained by the magnification 
peaks observed during a caustic crossing. If the light curve 
only samples up to some minimum physical distance from a 

9 In Einstein units, one can achieve the observed variability using either a 
large source moving rapidly or a small source moving slowly, with a degen- 
eracy of roughly r oc v. We always impose a prior on the physical velocity 
V oc v(M) 1 / 2 , so the physical source size r oc r(M}'/ 2 oc i>(M)'/ 2 oc v is es- 
sentially independent of (M) given a prior on the velocity. 
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FIG. 5. — The X-ray (top) and R-band optical (bottom) flux ratios between 
the A— B and B— C images along with the tracks across the microlensing pat- 
terns for images A (left) and B (right). The large circle shown on each pattern 
is the Einstein radius, while the small circles have the half-light radius of the 
optical disk and are shown at the positions corresponding to the epochs of 
the X-ray observations. The overall length of the line corresponds to one 
decade of motion. Darker colors represent logarithmically higher magnifica- 
tions with an overall magnification range from 1/30 to 30. This is a Case 
I example with fairly large differential offsets. It has a high stellar surface 
density (/* = 0.7), a large amount of smooth optical emission (/hop = 0.4), 
and the X-ray source is 14 times smaller than the optical. 

caustic crossing, then it will constrain sources sizes signif- 
icantly smaller than that distance poorly and the likelihood 
function will flatten for small source sizes. A logarithmic size 
prior then favors these small scales compared to a linear prior, 
leading to significant differences. Thus, our lower limits on 
the size of the X-ray emission are at a minimum prior depen- 
dent. More conservatively, the results could be interpreted as 
providing only an upper bound on the size of the X-ray emit- 
ting region. 



FIG. 6. — As in Fig. [5] This is a Case IT example, so the magnification 
offsets are small and the same for both the optical and X-ray data. It has a 
very low stellar surface density (/* = 0.1), a little smooth optical emission 
(/nop = 0.2), and the X-ray source is 32 times smaller than the optical. 

Fig EJshows how the sizes depend on the priors, the treat- 
ment of the magnification offsets and the fraction /nop of the 
optical emission that is unaffected by microlensing. The X- 
ray size is affected only by the choice of the size prior. The 
optical size is only affected by /nop- There are no significant 
differences between the results for the four magnification off- 
set cases. In order to produce the same optical variability with 
a larger fraction /nop> we must shrink the disk scale length 
R\. Mortonson et al. (2005) argue that the effects of mi- 
crolensing are largely determined by the half-light radius of 
the source, which is R\n = 2.44R\ in the limit that /nop = 0. 
As we increase /nop, the disk scale length shrinks roughly 
by the amount needed to keep the half light radius constant, 
with R x /2 = 5.2HR\ when /nop = 0.4. Note, however, that 
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FIG. 7.— AsinFig.[5] This is a Case F example, so the magnification 
offsets are small but differ between the optical and X-ray data. It has a low 
stellar surface density (/* = 0.3), a little smooth optical emission (/no/j = 
0.2), and the X-ray source is 28 times smaller than the optical. In the left 
panel we are at the edge of the pattern (although the Kochanek (2004) pattern 
creation method here produces periodic patterns that allow us to wrap the 
light curves across edges). 

the scaling for this particular model will break down when 
fno^ = 1/2. 

The larger values of /noju are favored, with likelihood ratios 
of 0.35, 0.40, 0.64, 0.92 and 1 .0 for / n0A1 = 0, 0. 1, 0.2, 0.3 and 
0.4 respectively. These differences are only marginally signif- 
icant, but they are in the sense of favoring (effectively) a flatter 
temperature profile. A flatter temperature profile can help to 
reconcile the differences between the larger microlensing and 
thin disk theory sizes as compared to the smaller flux sizes. 
Such flatter temperature profiles are generally consistent with 
the studies of the optical/infrared wavelength dependence of 
microlensing (Anguita et al. 2008, Eigenbrod et al. 2008, 



Poindexter et al. 2008, Mosquera et al. 2009, Bate et al. 
2009), but are not required. The one exception is Floyd et al. 
(2009), who find a limit requiring a steeper temperature pro- 
file. Some of this information is also present in the overall 
spectral energy distribution, and it is a long standing problem 
that the spectra of quasars do not match the predictions of thin 
disk theory (see Koratkar & Blaes 1999, Gaskell 2008). 

Whether increasing fno^ helps to resolve the size discrep- 
ancies depends on the physical model for the contamination. 
Line emission is reprocessed shorter wavelength emission, 
so as we increase /no^ we are a ls° reducing the fraction of 
the observed emission due to the disk and the flux size also 
shrinks as (1 -/no^) 1//2 - If, however, we view it as scattering 
fraction fno^ of the continuum emission on some large scale, 
then the flux size estimate is unchanged and the effect helps to 
reduce the discrepancy. Resolving the discrepancy with /no/j 
would require that most of the optical emission does not reach 
us directly from the accretion disk. 

While the source sizes show little dependence on the treat- 
ment of the magnification offsets, estimates of the amount of 
dark matter in the lens are sensitive to how strongly we con- 
strain the models to match the observed macro model flux 
ratios, as illustrated in Fig. @] By leaving the offsets rela- 
tively free, so as to conservatively estimate the source sizes, 
we have not optimized the calculation for probing dark matter. 
The Case I and II models, where we very loosely constrain 
the allowed magnification offsets, marginally favor models 
with /* ~ 0.3. The Case F and W models, where we only 
accept small offsets, favor the same dark matter dominated 
model more strongly. This range for /* is also that favored by 
the time delays measurements from Kozlowski et al. (2009). 
Note, however, that we are not in the "lagoon and caldera" 
regime for the microlensing patterns noted by Schechter & 
Wambsganss (2002) because of the very high magnifications 
(/i ~ 50 for image A at low k*/k rather than fi ~ 10). 

Figs. EE] show several examples of model light curves that 
fit the data reasonably well. These were also selected to have 
velocities consistent with masses of order 0.1-1.0M Q . The so- 
lutions are not unique, but they illustrate how simple changes 
in the source size dramatically alter the amplitude of the vari- 
ability. A common theme to the solutions is that the A and B 
images are generally required to lie in "active" regions of the 
patterns in order to produce such large changes in the X-ray 
fluxes. This means that we can expect the dramatic variability 
observed in this system to continue for an extended period of 
time (1-10 years). It is also interesting to note that significant 
changes in the optical fluxes are also likely. The larger op- 
tical source size washes out the effects of the closely spaced 
caustics that help to drive the X-ray variability. But the over- 
all changes between the high magnification ridges and the de- 
magnified valleys are still significant, and we should see over- 
all changes in the optical fluxes several times those observed 
to date. 

The implications of these results for theoretical models are 
mixed. The size of the disk is grossly similar to that expected 
from thin disk theory, and as we have summarized in Morgan 
et al. (2009), we find disk sizes that scale with black hole 
mass and optical wavelength roughly as expected (Eqn. [5]). 
We also find that the X-ray emission arises from significantly 
closer to the expected inner edge of the accretion disk than the 
optical emission, as we might expect for a hot corona. The 
optical size is broadly consistent with the expectations for an 
Eddington luminosity black hole with a mass estimated from 
the emission line widths (Eqn. [5]). But the size is inconsis- 
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tent with that expected for a thermally radiating disk with the 
observed magnitude (Eqn. |3). This is the discrepancy origi- 
nally noted by Pooley et al. (2006), which we explore more 
quantitatively in Morgan et al. (2009). 

Should we conclude that the thin disk model is wrong or 
simply that we have oversimplified the optical radiation trans- 
fer? We considered contamination by line emission or scat- 
tering of the optical photons, finding that this can modestly 
reduce the disk size for the range where up to 40% of the op- 
tical emission does not come directly from the disk. Our sim- 
ple emission model neglects the disk atmosphere and heating 
of the outer disk by radiation from the inner disk, all pro- 
cesses which would tend to make the optical emission region 
larger than the point where the disk has a temperature match- 
ing the photon wavelength without any change in the under- 
lying properties of the disk. Many of these effects are in- 
cluded in recent disk models such as Hubeny et al. (2001) or 
Li et al. (2005) for disk spectra. We examined face-on mod- 
els with M BH = 10 8 M Q , M = 0.09M Q yr" 1 and a BH spin of 
a = 0.998 using the Hubeny et al. (2001) models to compute 
our definition of the disk scale (R\, where IcT = hc/X) and 
the half-light radius (Rin)- The scale R\ is the most sensi- 
tive to the assumptions, with R\/r g = 41 for \ t, s = 0.814/im 
in our simplified disk model but equal to 36/34 for black 
body non-relativistic/relativistic disk models (BB NR/REL) 
and to 28/26 for non-LTE non-relativistic/relativistic models 
(NLTE NR/REL). The model dependence is much reduced if 
we compare with the half-light radii (R\/2/r g = 100 for the 
simple model, 114/117 for the NLTE NR/REL models, and 
99/102 for the BB NR/REL models). A flatter temperature 
profile than T oc /?" 3 / 4 would help, since at fixed total flux 
the half light radius increases. For example if T oc R~ 1 / 2 , the 
flux would be only 20% that of our standard profile for the 



same half-light radius. Indeed, such a flat temperature profile 
would also come much closer to matching the observed spec- 
tra of quasars (e.g. Koratkar & Blaes 1999, Gaskell 2008) and 
would be representative of models dominated by irradiation. 
In general, however, current microlensing results on tempera- 
ture profiles do not favor such flat profiles even if they gener- 
ally allow somewhat flatter profiles (e.g. Anguita et al. 2008, 
Eigenbrod et al. 2008, Poindexter et al. 2008, Mosquera et al. 
2009, Bate et al. 2009, Floyd et al. 2009). 

The key to disentangling these problems is to expand the 
measurements over a broad range of wavelengths, so that we 
can constrain the temperature profile of the disk, and over a 
broad range of black hole masses and accretion rates. For 
the particular case of RXJ1131, we have programs to con- 
tinue the X-ray monitoring of the system and to use HST to 
monitor the ultraviolet flux ratios of the images. Obtaining 
a robust lower limit to the size of the X-ray emitting region 
may require denser sampling of the X-ray light curve. Mea- 
suring the mid-infrared flux ratios of the images would be the 
most important step towards rigorously imposing the relative 
macro magnifications with no offsets, as these would only be 
affected by the macro model and any large scale gravitational 
perturbations (satellites). 
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TABLE 1 

MlCROLENSING MODEL PARAMETERS 



/. 


Image 


K 


7 




0.1 


A 


0.667 


0.359 


0.030 




B 


0.631 


0.325 


0.027 




C 


0.644 


0.306 


0.028 




D 


1.079 


0.493 


0.079 


0.2 


A 


0.618 


0.412 


0.067 




B 


0.581 


0.367 


0.060 




C 


0.595 


0.346 


0.062 




D 


1.041 


0.631 


0.159 


0.3 


A 


0.569 


0.465 


0.110 




B 


0.530 


0.410 


0.099 




C 


0.546 


0.387 


0.103 




D 


1.001 


0.635 


0.242 


0.4 


A 


0.519 


0.518 


0.162 




B 


0.480 


0.453 


0.146 




C 


0.496 


0.427 


0.153 




D 


0.964 


0.895 


0.329 


0.5 


A 


0.469 


0.572 


0.226 




B 


0.430 


0.497 


0.204 




C 


0.447 


0.469 


0.214 




D 


0.925 


1.018 


0.421 


0.6 


A 


0.419 


0.626 


0.305 




B 


0.379 


0.540 


0.278 




C 


0.397 


0.511 


0.290 




D 


0.890 


1.139 


0.520 


0.7 


A 


0.369 


0.679 


0.406 




B 


0.329 


0.584 


0.375 




C 


0.348 


0.552 


0.390 




D 


0.851 


1.288 


0.625 


0.8 


A 


0.318 


0.734 


0.541 




B 


0.278 


0.628 


(L507 




C 


0.297 


0.595 


0.524 




D 


0.816 


1.412 


0.740 


0.9 


A 


0.268 


0.787 


0.725 




B 


0.228 


0.671 


0.697 




C 


0.247 


0.637 


0.711 




D 


0.781 


1.530 


0.863 


1.0 


A 


0.217 


0.842 


1.000 




B 


0.178 


0.714 


1.000 




C 


0.196 


0.679 


1.000 




D 


0.740 


1.639 


1.000 



NOTE. — The macro models are parametrized by /„, the fraction of mass in the de Vaucouleurs component relative to a constant M/L model with /* = 1. 
The microlensing model parameters are the convergence k, shear 7 and the fraction of the convergence in stars k,/k. 



